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A complex scalar field has recently been suggested to bind galaxies and flatten the rotation curves 
of spirals. Its cosmological behavior is thoroughly investigated here. Such a held is shown to be 
a potential candidate for the cosmological dark matter that Hlls up a fraction flcdm ~ 0.3 of the 
Universe. However, problems arise when the limits from galactic dynamics and some cosmological 
constraints are taken simultaneously into account. A free complex held, associated to a very small 
mass m ~ 10“^® eV, has a correct cosmological behavior in the early Universe, but behaves today 
mostly as a real axion, with a problematic value of its conserved quantum number. On the other 
hand, an interacting held with quartic coupling A ~ 0.1 has a more realistic mass m ~ 1 eV 
and carries a quantum number close to the photon number density. Unlike a free held, it would 
be spinning today in the complex plane - like the so-called “spintessence”. Unfortunately, the 
cosmological evolution of such held in the early Universe is hardly compatible with constraints from 
nucleosynthesis and structure formation. 


I. INTRODUCTION 


The observations of the Cosmic Microwave Background (CMB) anisotropies [1], combined either with the determi¬ 
nation of the relation between the distance of luminosity and the redshift of supernovae SNela [2], or with the large 
scale structure (LSS) information from galaxy and cluster surveys [3], give independent evidence for a dark matter 
density in the range = 0.13 ± 0.05 [1], to be compared to a baryon density of Ub = 0.019 ± 0.002 as 

indicated by nucleosynthesis [4] and the relative heights of the first acoustic peaks in the CMB data. The nature 
of that component is still unresolved insofar. The favorite candidate for the non-baryonic dark matter is a weakly- 
interacting massive particle (WIMP). The so-called neutralino naturally arises in the framework of supersymmetric 
theories. Depending on the numerous parameters of the model, its relic abundance Dcdm falls in the ballpark of the 
measured value. New experimental techniques have been developed in the past decade to detect these evading species. 
However, detailed numerical simulations have recently pointed to a few problems related to the extreme weakness 
with which that form of matter interacts. Neutralinos tend naturally to collapse in numerous and highly packed 
clumps [5] that are not seen - see however [6]. The halo of the Milky Way should contain ~ half a thousand satellites 
with mass in excess of 10® Mq while a dozen only of dwarf-spheroidals are seen. The clumps would also heat and 
eventually shred the galactic ridge. More generally, this process would lead to the destruction of the disks of spirals. 
A neutralino cusp would form at their centers. This is not supported by the rotation curves of low-surface-brightness 
galaxies that indicate on the contrary the presence of a core with constant density. Finally, two-body interactions 
with halo neutralinos and its associated dynamical friction would rapidly slow down the otherwise observed spinning 
bars at the center of spirals like M31. 

New candidates for the astronomical dark matter are therefore under scrutiny such as warm dark matter [7], particles 
with self interactions [8], or non-thermally produced WIMPs [9]. An exciting possibility would be to have a common 
explanation for both the dark energy Da and the dark matter Dcdm components of the Universe. Before trying to 
reach such an ambitious goal^, one could explore the relevance of scalar helds to the cosmological and galactic dark 
matter puzzles, as was done for dark energy with the so-called “quintessence” models [11]. The archetypal example 
of quintessence is a neutral scalar field ip with Lagrangian density 
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( 1 ) 


^ = ^9^''d^ipdu(p - V {(fi) . 

Should the field be homogeneous, its cosmological energy density would be expressed as 

P^ = Tl = ^ +Viv’) , ( 2 ) 

whereas the pressure would obtain from = —gij P so that 

= Y ~ 

If the kinetic term ip'^/2 is small with respect to the contribution from the potential V {(p), the equation of state can 
match the condition for driving accelerated expansion in the Universe, loq = P^p/Pip < —1/3. Instead, in order to 
behave as dark matter today, the field should be pressureless: \Pp\ <C Pp. So, the kinetic and potential energies should 
cancel out in Eq. (3), a condition automatically fulfilled by a quickly oscillating scalar field averaged over one period 
of oscillation. This well-known setup is that of the cosmological axion. It requires a quadratic scalar potential, so 
that the kinetic and potential energies both redshift as oc a~^ with the Universe expansion and cancel out at any 
time during the field dominated stage, which is then equivalent to the usual matter dominated one. 

Axions - or more generally, bosonic dark matter - were revived recently due to the undergoing CDM crisis. For 
instance, it was noticed in [12] that stucture formation on small scales can be forbidden by quantum mechanics, for 
wavelengths smaller than the Compton wavelength - i.e., the minimal spreading of an individual boson wave function. 
The latter matches the scale of galactic substructures for an ultra-light mass of order m ~ 10“^^ eV. Alternatively, 
one may introduce a self-coupling term [13,14]. As we have seen, the existence of a matter-like dominated stage 
requires that the contribution of non-quadratic terms to the potential energy remains subdominant. Nevertheless, 
a self-coupling would modify the field behavior in the early Universe, as well as its clustering properties today in 
regions where the field is overdense - exactly like for boson stars, which are crucially affected by the presence of a 
self-coupling [15]. The self-coupling is also relevant for the issue of field clumps stability, and can explain why dwarf 
and low-surface-brightness galaxies have cores with finite density [16]. 

A remarkable feature with bosonic dark matter is the possibility to form Bose condensates, i.e., large domains 
where the field is coherent in phase and is in equilibrium inside its own gravitational potential - like boson stars - 
or in that of an external baryonic matter distribution. This opens the possibility to have a very simple and elegant 
model for galactic halos, in which rotation curves would follow from simple equations - essentially the Klein-Gordon 
wave equation, which governs classical scalar fields as well as Bose condensates. This situation strongly differs from 
the more conventional picture of a gas of individual particles - fermions or heavy bosons - for which gravitational 
clustering does not lead to universal density profiles and where the shape of galactic halos can only be studied through 
technically difficult A^-body simulations. 

The formation and stability of such condensates is a complicated issue - see e.g. [19-21] - even when the field is 
complex and has a global charge - not to be understood as an electric charge, but as a conserved number of quanta 
like the baryon or lepton number. For instance, a large condensate can be unstable under fragmentation into smaller 
clumps. For a real scalar field with no conserved charge, the issue of stability is even more subtle since the field can 
self-annihilate, especially when the condensate core density exceeds a critical value [19]. This property can improve 
the agreement with observations [16], since the coupling constant will tune the upper limit on the density of dark 
matter cusps at the center of galaxies. However, such a positive feature is far from excluding models with a conserved 
charge. In fact, the issue of Bose condensation on galactic scales - in an expanding Universe - has never been studied 
in details. The result would depend very much on the scalar potential, and it is difficult to guess what would be the 
maximal core density today. 

In this work, we will focus on the scenario with a conserved charge, and assume that dark matter consists in a 
complex scalar field with a quasi-homogeneous density in the early Universe, producing later galactic halos through 
Bose-condensation. The Lagrangian density reads 

C = - V {cfy . (4) 

Throughout this analysis, we will focus on the potential 

V = + A . (5) 

As a prologue to the study of density perturbations, we will follow the evolution of the homogeneous cosmological 
background of this field, taking into account the constraints on the scalar potential coming from galactic halos. In a 
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previous work [17], we did a detailed comparison of such halos with universal galactic rotation curves [23], for a massive 
complex scalar field. We recall the salient features of this analysis in section II and discuss how the corresponding 
constraint on the mass m is modified when a self-interaction coupling A is introduced. In the subsequent section III, 
we study how this scenario can be implemented cosmologically. Neglecting the spatial variations of (j), we are led to 
the cosmological density 

+ V {4>) (6) 

and pressure 

. (7) 

Beside the issue of charge conservation, the case for a complex scalar field is somewhat richer than that of a real 
(neutral) scalar field. In one limit, the complex field can behave as an effective real one, similar to the usual axion. 
On the other hand, it can be spinning in the complex plane with slowly-varying modulus, like in the so-called 

spintessence [18] scenario. This depends on the dominant term in the kinetic energy, which can be either radial and 

oscillating, or orthoradial and slowly varying. As a result, during the field dominated era, the spintessence would 
have a continuously vanishing pressure, while the axion pressure would oscillate between and —p^. Finally, in 
section IV, we discuss our results and feature the problems that plague the axion-spintessence dark matter scenario. 
We will suggest some future directions worth being investigated. 


II. GALACTIC HALOS 


We are interested in galactic halos consisting in self-gravitating scalar field configurations - which can be seen as 
Bose-Einstein condensates spanning over very large scales. Since the typical velocities observed in galaxies do not 
exceed a few hundreds of km.s“^, it is enough to study the quasi-Newtonian limit where the deviations from the 
Minkowski metric = diag(l, — I, — I, — I) are accounted for by the vanishingly small perturbation h^n. Inside 
galaxies, the latter is of order the gravitational potential 

, ( 8 ) 


where Vc is the rotation velocity - in the case of spirals - and where Vq is the escape velocity from the system. 
In the harmonic coordinate gauge where it satisfies the condition 

= 0 , (9) 

the perturbation is related to the source tensor 

( 10 ) 


through 


= - I&ttQ Sfjtv 


( 11 ) 


The Poisson equation reads like 


A$ = SttQ ^oo 


( 12 ) 


where $ = /ioo/2 is the Newtonian potential. For pressureless matter, 2 Sqq = Too = P- On the other hand, for the 
complex scalar field, the gravitational potential is sourced by the effective mass density 


eff 

— c 

— = i’oo - 


{2<^t(^ - V{<P)} 


(13) 


which is a priori different from its cosmological counterpart (6). So, inside a galactic halo, the gravitational potential 
is given by 


A$ = 47r^ {p’f + pb) 


(14) 
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where pb is the distribution of baryonic matter forming the various galactic components - stellar disk, bar, bulge... 
In first approximation, the galaxy can be seen as spherically symmetric. In that case, one shows [22] that all stable 
field configurations must be in the form 


0(r, t) 


Jojt 

V2 


(15) 


where the amplitude a depends only on the radius r. Then, the effective field density reads like 


Ptf, 


o 2 2 z 

Zuj a — m a 


2^2 



(16) 


The radial distribution of the field a{r) and the gravitational potential <i)(r) are given by a system of two coupled 
equations: the Poisson equation (14) and the Klein-Gordon equation. The latter may be expressed as 


e S a 


I 2 —2u 2 \ 3 n 

+ uj e a — m a — Act = 0 


(17) 


in the isotropic metric where 

dr^ = dt^ - {dr^ + r'^dO^ + sin^ 0 V} . (18) 

The prime denotes the derivative with respect to the radius r. The Newtonian approximation corresponds to 


u ~ — u ~ 4) . (19) 

Relation (17) simplifies into 

Act + (1 — 44>) a — (1 — 24>) {n? a + Act^) = 0 . (20) 

For each value of the parameters (m, A, cu) and a given baryon distribution, these equations form an eigenvalue 
problem with discrete solutions, labelled either by the central value ctq > 0 or by the number of nodes n in which 
a(r) = 0. The lowest energy state - which is not identically null due to charge conservation - has n = 0. The 
self-consistency of the Newtonian limit requires |4>| <C 1. Such solutions exist only for 

0 < (m^ — uj^) . (21) 


A. Free field 


In [17], we solved these equations for A = 0. We found that halos consisting in the fundamental configuration of 
a free scalar held ht perfectly well the universal rotation curves of low-luminosity spiral galaxies [23]. These data 
has three advantages for our purpose: the robustness of the points and error bars - obtained by averaging over many 
galaxies -, the good determination of the baryon distribution - solely a stellar disk with exponential luminosity prohle 
- and the low baryon contribution which justihes the approximation of spherical symmetry. 

With a quadratic potential, the size of the halo is given by 


I 


jMp 

h 

/ O-Q 

me 


( 22 ) 


where we neglected the dependence on the baryon density. If the central held value ctq is signihcantly smaller than 
the Planck mass, the coherence length of the condensate exceeds the Compton wavelength of an individual particle - 
^Compton = ft/{me) - but it is clear that only an ultra-light scalar held can condensate on distances of order 10 kpc. 
The typical orbiting velocity in such a halo is given by v/c ~ ao/Mp. Therefore, requiring v ~ 100 km.s“^ and 
I ~ 10 kpc hxes ctq around 10”® Mp and m around 10”^® eV, as conhrmed by a detailed htting to the data [17]. 

Since the distribution of such halos only depends on the free parameters ctq and m - where we impose a unique 
value of m for all galaxies - we believe that their success in reproducing universal rotation curves is a signiheant 
argument in favor of this model. On the other hand, the existence of such a low mass, even if not strictly forbidden 
by fundamental principles, is very unlikely due to unavoidable radiative corrections. This could motivate a systematic 
investigation of other potentials for the scalar held. The next level of complexity would consist in adding a quartic 
self-coupling. 


4 





B. Quartic self—coupling A 


As is well known for boson stars - which are exactly identical to our halos in the absence of a baryon component 
- the inclusion of a quartic term drastically modifies the mass m and the extension of the condensate, even when it 
contributes to a negligeable fraction of the central energy density [15,24]. This is so because Acr^ can be very small 
with respect to ~ and yet comparable to the difference (m^ — that appears in the Klein-Gordon equation. 
In the limit where A = A/(47r0m^) ^ 1 and in the absence of a baryon population, we can even give an exact 
analytic solution for the field and for the orbiting velocity of test particles: 



with the requirement that 

A-i « ^ « A-i/2 . (25) 

Because a^/Mp ~ (r;/c)/-\/A, the second inequality follows from the Newtonian self-consistency condition |<I)| ~ ^ 

1. The first inequality translates into A ^ (c/u) . It implies that all the field spatial derivatives can be neglected 
in Eq. (20) and sets the maximal radius up to which the analytic solution is valid. This maximal radius is at most 
equal to half a period so that r < y^A/2 (7r/m). That bound is almost saturated for A —> oo. Note that w does not 
appear in the analytic solutions because it is only relevant at larger radii. However, the Newtonian self-consistency 
condition imposes that I — uP' jw? ^ I. 

The field behavior (23) may be readily recovered by neglecting the spatial derivatives of cr in the Klein-Gordon 
equation (20) so that 


m 


02 


"'W = - (-8 - ‘ 

where B(r) =1 + 2 $(r) and O = w/m. In the Newtonian limit, the pressure (7) reads like 


P^ = C - 


4A 1 B 


- 1 


whilst the effective mass density (16) of the Bose-condensate is 


Ptf> 


- m ^ -I 


Both are related through the Lane-Emden polytropic equation of state 

B<j> — B P(j) 1 


(26) 


(27) 


(28) 


(29) 


with K = X /(4m^) while the polytropic index is n = I. For such a value, the gravitating system - in hydrostatic 
equilibrium - is shown to have a constant core radius a where 


a2 = 


I A 


8ttQ mX 

The field and density profiles are functions of the reduced radius z = rja 

P' 4 ,{r) _ a^(r) _ sin( 2 :) 

P0(O) erg 


(30) 


(31) 
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The most striking feature in the large A limit is as follows: although the quartic term remains subdominant in the 
energy density - Eq. (25) implies that \cj'^ <C - the typical size of the system is very different from the free field 

case since now it reads like 


. (32) 

As the central field value does not appear in this expression, different halo sizes would just result from different baryon 
contributions to the density, which bring corrections to (32). The central field value ao still determines the rotation 
curve amplitude. In the large A limit and in the absence of baryons, the maximal rotation speed is given exactly by 

2 2 

^ = 2.13 TT A ^ at r = 1.94A^/^/m . (33) 

In this paper, we will not extend our previous comparison with galactic rotation curves [17] to the case of a quartic 
self-coupling, leaving this for a future publication. We will just use relations (32) and (33) in order to find a rough 
order of magnitude estimate for A/m^ which has the best chance to provide a good ht to the data. In the next section, 
this constraint will be plugged into a cosmological dark matter scenario. By requiring that the rotation velocity peaks 
around 200 km.s“^ at a typical radius of 10 kpc, we find 

m ~ eV and (Tq ~ A“^/^ 10“^ Mp . (34) 

Taking for instance A in the range [1, 10“"^], we obtain a mass of order 0.1 to 1 eV, z.e., a few orders of magnitude 
larger than the expected neutrino masses. This gives A ~ 10^^ and (Jq ~ 0.1 eV. 

III. COSMOLOGICAL DARK MATTER 


A. Matter—like behavior 


A homogeneous complex scalar field with a quadratic potential is a perfect candidate for pressureless cold dark 
matter. Let us write the Klein-Gordon and Friedman equations for the field configuration 


m = . 


This amounts to 


cPcr 3 da da n , o o 

H-rr^ + + Acr'^ - uj^a = 0 

a dt dt 

dhj 3 da da 

—a H - —Loa + 2u!— = 0 

dt a dt dt 


da A 2 2 I 2 2 I 4 

-) +u;a +ma +-a 


(35) 

(36) 

(37) 

(38) 


where prad is the usual density of relativistic photons and neutrinos and uj = dO/dt. The second equation implies the 
conservation of the charge per comoving volume Q = coa'^a^. Therefore, we can rewrite the first equation as 


d^cr 3 da da , , o 


It is most convenient to use the dimensionless time and field variables: 

t = mt , H = H/m, u 

Then, the charge conservation implies lou^ = m and the Klein-Gordon equation reads like 


/m 3/2 

r "• 


u -\- {1 — u ‘^ + 


XQ 


H 


-H‘ 


u = 0 


(39) 


(40) 


(41) 
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where the dot denotes a derivative with respect to i. When A 7 ^ 0, we see that the self-coupling term always dominates 
the mass term in the past, when a —> 0. However, we will first study the opposite case when the quartic term is zero 
or subdominant. If Eq. (41) is to be applied today, or in the late stage of evolution of the Universe, we can also neglect 
the contribution from the Universe expansion: today, one has Hq ~ many orders of magnitude below the 

values of m discussed in section 11. So, Eq. (41) reduces to 

u + u — u~^ = 0 , (42) 

and describes some periodic oscillations in a static potential V(u) = {v?' + u“^)/2, with a minimum at mq = 1 
corresponding to uj = m. The field density reads 


= mQ EuQ ^ , (43) 

where £!„ > 1 stands for the conserved energy (v?+u^ +u~'^)/2. We conclude that in the late evolution of the Universe, 
as soon as m becomes bigger than H and A^/^cr, the homogeneous scalar field behaves exactly as a cosmological 
background of dark matter. 

Let us provide further intuition on the physical meaning of the oscillations for the variable u. If Eu — 1 and u is 
almost stable around one, then the modulus slowly decreases as tr oc while the phase velocity is constant with 

Lu = m. The equation of state is that of pressureless matter: 

P 

w=^ -J < 1 . (44) 


Such a field, following a spiral trajectory, is a particular case of what was recently called spintessence [18]. On the 
other hand, when ^ 1 , u and uj strongly oscillates, but the trajectory of is a fixed ellipse as can be seen 

from the exact analytic solution to Eq. (42): 


{(u2-l)^/%in (2t + a) +£:„} 


1/2 


j-TT/m 

0 (t = tt) — 0(0) = / uj{t)dt = 

d 0 


-2 


dt = 


(45) 


The integral of the phase over one period of oscillation obtains from Gradshteyn and Ryzhik, 3.661.4. So, the field 4> 
follows an ellipse which axes decrease as In the large Eu limit, the ellipse is squeezed and reduces asymptotically 

to a line. Then, the complex field is essentially similar to a real oscillating field, like in usual axionic dark matter 
models. The pressure does not vanish identically, but quickly oscillates between and —P 0 , with zero average over 
one period At = 7 r/m. This period is much smaller than E[~^: on cosmological scales, the field still has the same 
effect as pressureless matter. 


B. radiation domination: the A = 0 case 


In the A = 0 case, we know that the mass should be of order m ~ 10“^^ eV, ten orders of magnitude above the 
present value of the Hubble parameter. It is easy to show that H will start to dominate over the mass at a redshift 
of order 6 x 10®. Therefore, the previous analytic solutions (43) and (45) only apply to the end of the radiation 
dominated stage. During radiation domination, the Klein-Gordon equation (41) simplifies into 


u + u + 



= 0 . 


(46) 


When H ^ m - i.e., El ^ 1 - the term u can be neglected in the above equation. We then obtain a simple non-linear 
equation for the variable v = 


- -T = 0 , 


(47) 


where the prime denotes derivation with respect to conformal time dt = a dr. There is an exact solution: 


C.+ 

m a \ a 


a2 / 


1/2 


(48) 
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Only two of the three constants (Ci, Cz) are independent since 


O 1 C 3 


4 


{a?Hf 


(49) 


Remember that a^H is also constant during radiation domination. The free parameters (Ci, C 2 ) can be conveniently 
expressed as a function of the initial conditions for the field at some initial time 


Cl 

C 2 

These results show that a quickly 
reads like: 


2 2(Ji{du / dt)i [du/dt)1 

= CT,- H- 77 -1- 77t -h 


— 2cLi 


m ' Hf Hf ^ 
iTi{da/dt)i (da/dt)1 wfcr. 






H? 


(50) 

(51) 


stabilizes at a value ^/Ui. During radiation domination, the field energy density 


C 3 , m 2 ^ ^ 


a2 


(52) 


Interestingly, we see that after some time the energy can be dominated by the contribution of the mass term, even if 
the latter can be neglected in the Klein-Gordon equation. 

We now understand what the generic evolution looks like, starting from an early time - say, for instance, at the 
end of reheating - at which the field density is smaller than the radiation density: first, decays as a“®; then, it 
stabilizes around the value rr?Cij2, and remains constant as long as H > m] finally, when H < m, the density decays 
as that of pressureless matter, takes over the radiation density, and drives the matter-like dominated stage. Such a 
scenario requires a single constraint: the constant value of the field density during radiation domination should be 
matched with the matter density extrapolated from today back to the time at which H = m: 

^ Cl-7x10-4 m2 . (53) 


The scenario is enough constrained to reach another conclusion: when the field density decays as matter, its dynamics 
is that of an effective oscillating real scalar field and not that of spintessence. Indeed, as soon as the density becomes 
constant during radiation domination, one has 


(a2M)2 C3 m2 
- 2 - 


91 

n 6 


< m^C? 


(54) 


On the other hand, a crude matching between the two expressions for the energy density (43) and (52) at the time 
when H = m gives: 


mQ Ey^a ^ — m^ Ci 




Et 


m^Cl 


(55) 


By combining Eq. (54) and (55), we find that i?2 ^ 1, which is the condition for the field trajectory in the complex 
plane to be a squeezed ellipse of constant phase. 

So far, we have assumed that the radiation density was initially dominant, in order to use the analytic solutions 
(48) and (52). Instead, if the field dominates initially - a situation that could be allowed before the time of Big Bang 
Nucleosynthesis (BBN) - the field density will nevertheless decay as a”®. Indeed, if > prad, then necessarily2 
P 0 ^ vn?a‘^. We conclude that = p^, and energy conservation implies that ptj, oc a”®. As soon as the radiation 
density takes over, the previously described analytic solutions do apply. 


2proof: even if a was as large as Mp, would be of order 10 which corresponds to the radiation density at 

redshift 10®, well after BBN. 













FIG. 1. Evolution of the field and matter densities, in the scenario with A = 0 and m ~ 10“^® eV. The evolution starts at a 
redshift z = 10^° - but we could have equally well started earlier. The initial value of the field density can be chosen arbitrarily, 
below or above that of radiation density. After decaying as a”®, the field density reaches a plateau which amplitude has been 
fixed according to Eq. (53). This condition ensures a correct value of the density today: = flcdmPc- At 2 = 6.3 x 10®, when 

H — m, the density starts to decrease as o“®, like for presureless matter, and takes over radiation at 2 ~ 3 x 10®. 


Q. 





FIG. 2. Two possible evolutions of the field real and imaginary parts, in the scenario with A = 0 and m ~ 10“^® eV. On 
the left panel, the initial held modulus has been chosen below the equilibrium value \/Ci as determined in Eq. (53). The 
opposite situation is shown on the right panel. In both cases, the held density decays like in Fig. 1. At redshift 2 = 6.3 x 10®, 
when H = m, the held starts to oscillate, but its density decays smoothly as that of presureless matter. Because the real and 
imaginary parts are exactly in phase, the held is equivalent to a single real scalar held. 


C. radiation domination: A 7 ^ 0 

If the field has got a quartic self-coupling, Atr'^ must be negligible with respect to in the late Universe in order 

to drive a matter-like dominated stage with oc a“®. However, a quartic self-coupling is likely to be cosmologically 
relevant at early times, whenever the field modulus a well exceeds In that case, the equation for v - see 

relation (47) - reads like 


XQ 


m"’ 



(56) 
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During radiation domination, a" = 0 and u is a periodic - elliptic - function, describing non-harmonic oscillations in 
the potential V{v) = \Qv‘^/(4m^) + 1/(2 z;^), with a minimum at vq = The period of the oscillations 

- expressed in conformal time - is of order . So, a performs damped oscillations with a constant period 

Aa ~ with respect to the scale factor. If we furthermore define the conserved pseudo-energy energy of 

V hy Ey = v'^/2 + V{v), we can express the field density as 

f Ha , H'^a^ 2 I -4 /rn\ 

= mQ < Ey - vv + -—> a ^ . (57) 

( TO 2 TO^ J 


Remember that Ha decays as a ^ during radiation domination: a priori, at early times, the field density performs 
damped oscillations like the field modulus while at late times it decays smoothly, as for radiation 

p^ = mQEya~^ . (58) 

The transition between both behaviors takes place when H^a^v^ ~ nri^Ey, where v is evaluated at the maximum of 
one oscillation: Vmax ~ (4m^Ey/AQ)^‘^^. Inserting this condition in Eq. (57), we find that the transition between 
damped oscillations and smooth decay occurs when p^ ~ H'^/X. In practice, this implies that the oscillatory behavior 
is generally irrelevant unless A is fine-tuned to extremely small values - for the ordinary radiation component, the 
condition prad ^ H^ is already realized at the end of inflation. So, Eq. (58) applies even in the early Universe. 

Later on, the transition between radiation-like and matter-like behaviors will be effective when the maximal value 
of a during one oscillation, computed from Eq.(56), will be comparable to X~^^^m. This translates into 


^max — 



1 / 4mQEy A TO 

a V ^ / VA 


( 4XQEy V'^ 

\ TO^ y 



(eV)4 . (59) 


In the last equality we used the constraint from the size of galactic halos. Since, on the other hand, peq ~ 0.55 (eV)^, 
the transition to matter-like behavior occurs slightly before equality. This means that at earlier times, when the field 
behaves like radiation, its density should be fine-tuned in order to be close to the radiation density. 

This scenario is illustrated on Eig. 3. In order to obtain the correct value of the field and radiation densities today, 
P 0 must be adjusted to 0.6 Prad- This would correspond to an effective number of extra neutrino species of AiVeff = 5 
that is not even allowed by BBN. 



a/ao = (1+z)''' 


FIG. 3. Evolution of the field and matter densities, in the scenario with m ~ A^^^ eV and arbitrary A. The evolution also 
starts at a redshift z = 10^°. The field decays first as dark radiation, and then as dark matter. The transition between these 
two regimes is constrained by Eq. (59) to take place immediately before matter-radiation equality. In other words, in the early 
Universe, the field density must be very close to that of radiation. The simulation gives p^, = 0.6 prad in order to obtain a 
correct value of p^ = pcdm today (such that flcdm ~ 0.13). 
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IV. DISCUSSION 


In section II, we recalled why a complex scalar field is an attractive candidate for dark matter in galactic halos; it 
can provide a rather powerful explanation to the observed rotation curves, because the radial distribution of a scalar 
field after Bose condensation is constrained by a wave equation (the Klein-Gordon equation); this form of dark matter 
would be very smooth, and unlike a gas of individual particles, it would not slow down the relative motion of the 
baryonic components through dynamical friction. 

The results of section III show that the complex scalar field can also play the role of a cosmological dark matter 
background, since after various possible behaviors in the early Universe, which depend on the exact shape of the 
potential, the field will decay like ordinary pressureless matter as soon as the mass term dominates the potential and 
the mass is bigger than the Hubble parameter. 

However, when the parameters of the scalar potential are estimated from the size and mass of galactic halos, and 
plugged into the cosmological evolution, some tension appears both for the free-held model and for the one with a 
quartic self-coupling. Indeed, the redshift at which the held starts to decay like matter, cx a“^, is completely hxed 
by these parameter values. 

In the case with A = 0, this transition occurs at a redshift close to z = 6 x 10®. The comoving wavenumber of 
perturbations entering into the Hubble radius at that time is fe = 3 Mpc“^. As shown in Ref. [12], this is of the same 
order of magnitude as the Jeans length for the free held. So, larger wavelengths - in particular, those probed by the 
spectrum of CMB anisotropies and by the linear matter power spectrum - are expected to undergo the same evolution 
as in an ordinary CDM model. On smaller wavelength, one would still need to prove with numerical simulations that 
Bose condensation can occur on the scale of galactic halos - and eventually also on slightly bigger and smaller scales. 
So far, this model still sounds attractive, apart maybe from the ultra-light mass required (m ~ 10“^®eV). 

In the case with A 7 ^ 0, the transition between radiation-like and matter-like behaviors happens even later, just 
before equality; therefore, the held density during radiation domination has to be very close to that of photons and 
neutrinos. Our simulation reveals that ^cdmh'^ = 0-13 today obtains from p^/prad — 0-6 during radiation domination; 
in other words, the held contributes to the number of relativistic degrees of freedom as an effective neutrino number 
AA"eff — 5, in contradiction with the BBN constraint |AA"eff| < 1- Moreover, the perturbations which enter into the 
Hubble radius before the transition - which would be probed today by the spectrum of CMB anisotropy and the linear 
matter power spectrum, since their comobile wavenumber is bigger than k ~ 10 “^ Mpc“^ - would be suppressed with 
respect to CDM perturbations (essentially, like for hot dark matter). These problems were previously noticed by 
Peebles [13] in the case of a real scalar held. As a possible solution, Peebles suggests a polynomial self-coupling with 
a power slightly smaller than four. In fact, in order to save this scenario, it would be enough to increase the redshift 
of the transition by a factor ten or so. We conclude that the model with a quartic self-coupling is marginally excluded 
by cosmological constraints, but that any small deviations from the simple framework studied here are worth being 
investigated. 

Our analysis is still incomplete. We said that the conservation of the charge Jo = — 4>4’*) was crucial for 

the stability of the condensates. Charge conservation also provides further constraints between, on the one hand, the 
cosmological background of the held in the early Universe, and on the other hand, the current distribution of dark 
matter in the form of held condensates. Also, it can give some hints on possible mechanisms for the generation of 
the held density in the early Universe. In terms of quanta, the conservation of charge implies a constant number of 
bosons minus anti-bosons inside a comoving volume: 

Q= n^)a^ = cte , (60) 

where (n^, n^) are the number density of bosons and of their CP conjugates. We want to evaluate this number today 
in the form of galactic dark matter. A priori, the rotating phase inside each condensate is hxed up to an arbitrary 
sign: 6{t) ~ ±mt. So, some halos can be made of bosons, and some others of anti-bosons. This may occur, for 
instance, if the scalar held was populated during a phase transition: after the transition, the initial held distribution 
could be very homogeneous, but with an arbitrary phase distribution, leading to domains with positive and negative 
values of w = dO/dt. A priori, such a disordered initial conhguration - with no phase coherence - may not exclude 
the subsequent formation of Bose condensates. In this case, the mean number density of bosons minus anti-bosons 
today can be arbitrarily close to zero, and the conservation of charge does not give any new constraint. Note however 
that two halos of opposite charge could annihilate and radiate out a massive amount of energy. Although this issue 
would deserve a more careful study, it is probably in conhict with observations. 

Let us consider now the case where all galaxies carry a charge with same sign. This can occur if the held underwent 
inhation - or was coupled to a held that underwent inhation - in such way to be very homogeneous in the early 
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Universe, as assumed in section III. Then, at the beginning of the matter-like regime, the phase would be coherent 
even on super-horizon scales, and all condensates would form with the same rotating phase. We can estimate the 
charge in galaxies, Qgai, by multiplying the typical charge of a single halo by the number density of halos. This charge 
must be smaller than or equal to the total charge of the cosmological homogeneous background, Qtot = 

In the A = 0 case, the charge inside one halo is given by [22,17]: 


N 



Ml 

m? 


10^® . 


(61) 


Under the very crude assumption that the Universe contains in average one halo per volume V = 1 Mpc®, we find a 
mean density today 

= Ogai Oq 3 ~ N/V ~ lO^^cm-® . (62) 

This number is extremely large, 10®^ times bigger than the present number density of photons. This is a strong 
indication that the model is not realistic: it would be very difficult to generate such a huge charge in the early 
Universe. In fact we can even completely exclude the model by comparing Qgai with the total charge of the cosmological 
background. We saw in section III, Eqs.(54, 55), that the existence of a plateau for during radiation domination 
imposes today a field dynamics close to that of an oscillating real axion, rather than spintessence. For spintessence, 
the number density = Qtot is equal to p^/m as can be seen from Eq.(43) with = 1. For an oscillating 
axion, most of the kinetic energy is in the radial direction and <IC p^/m since the term in relation (43) is now 
much larger than 1. Because we must be in the latter case at least during matter domination, we find the following 
upper bound on the total charge today: 

= Qtot aQ^ < - 10®^cm“® . (63) 

We hnd Qtot ^ Qgai, which is impossible. We conclude that the present dark matter scenario, based on a complex 
free scalar field forming galactic halos after Bose condensation, is not consistent - at least when the field is assumed 
to be homogeneous in the early Universe, and today all condensates carry a charge with the same sign. It would 
certainly be interesting to investigate the opposite scenario, with a homogeneous initial density but random phases, 
with the drawback that halos and “anti-halos” may annihilate. 

As we have seen, the case A 0 is already marginally excluded by cosmological constraints, but it is worth 
calculating the various charges also for this model. An individual halo has got a charge - see relation (34) - 


\m 


o-Q 



Under the assumption that all halos have a positive charge, one finds 


^<p — Qgal ^ 


-3 

0 


N/V - 30 cm"^ 



(64) 


(65) 


We are led to two intriguing coincidences. First, with A of order one - and therefore m of order 1 eV - the field 
number density is of the same order of magnitude as that of photons for which Uj = 400 cm“^. Second, provided 
that the field behaves cosmologically as spintessence, with p^ m ricf,, we can calculate the total charge and find 
Qtot ~ Qgai: the cosmological and astrophysical charges are consistent. Therefore, the scenario requires a mechanism 
in the early Universe that would fine-tune both and to some values very close to and pj. 

As we already said, this scenario with a quadratic coupling is plagued by inherent difficulties to produce small-scale 
perturbations, and by an inconsistency with the number of relativistic degrees of freedom indicated by nucleosynthesis. 
However, one should retain two positive features. First, any small modification of the scenario that would shift by a 
factor ten the redshift of the transition between radiation-like and matter-like behaviors would evade these difficulties. 
Second, the initial conditions for the field number density and energy density should be surprisingly close to those of 
photons. 


Throughout this discussion, we tried to give some arguments both in favor and against the two models considered 
here, based on two different scalar potentials. In their present form, none of these models can survive. However, we 
believe that one should retain the many positive indications discussed before as an encouragement for investigating 
other variants of scalar field dark matter. The fact that the two scalar potentials lead to very different conclusions 
already shows how rich and unpredictable is the field. 
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